Size scaling of strength in thin film delamination 
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Abstract. 

We investigate by numerical simulation the system size dependence of the shear 
I delamination strength of thin elastic films. The films are connected to a rigid substrate 

by a disordered interface containing a pre-existing crack. The size dependence of the 
, strength of this system is found to depend crucially on the crack shape. For circular 

P3 ' cracks, we observe a crossover between a size-independent regime at large crack radii 

O . which is controlled by propagation of the pre-existing crack, and a size-dependent 

regime at small radii which is dominated by nucleation of new cracks in other locations. 
For cracks of finite width that span the system transversally, we observe for all values 
of the crack length a logarithmic system size dependence of the failure stress. The 
' results are interpreted in terms of extreme value statistics. 
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1. Introduction 

The strength of disordered materials depends on multiple characteristic length scales, 
ranging from the macroscopic size of the sample over the sizes of pre-existing cracks down 
^ ■ to characteristic length scales of the disordered microstructure. Sample size dependence 
' is commonly related to the concept of a 'weakest link': Larger samples on average 
contain more weak parts and thus, provided failure is governed by the weakest part 
of the microstructure (e.g., the largest crack), they are likely to fail at lower stress. 
This argument puts the size dependence of failure firmly into the realm of extreme 
value statistics. Samples containing a single large crack, on the other hand, fall into 
an entirely different category: As failure is governed by the supercritical propagation of 
this crack, the behavior is expected to depend on the crack size but - at least in the 
large sample limit where boundary effects are negligible - not on the size of the sample. 

In disordered quasi-brittle materials, the above arguments must be used with 
caution. In such materials, large flaws can emerge from the collective dynamics of 
interacting microcracks. Thus, the concepts of 'weakest link' or 'critical flaw' are far 
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from straightforward - rather than being fixed pre-existing entities, weak zones and the 
associated length scales may need to be envisaged as emergent features of the complex 
dynamics of interacting defects. At the same time, the properties of pre-existing cracks 
are strongly influenced by their interactions with the disordered microstructure, which 
may give rise to crack front roughening [H [2] , the formation of self-affine fracture surfaces 
[3] , |1] and the correlated nucleation of microcracks in a process zone ahead of the crack 
tip [6]. 

The fracture behavior of structurally disordered materials has been studied in 
numerous investigations (for overview, see [5]). The size dependence of failure stresses 
in pre-cracked samples, which is in the focus of the present study, has been investigated 
for two-dimensional disordered sheets by Alava and co-workers f6j. In that study, the 
crack was one-dimensional and contained within the plane of the sheet. In the present 
study, we also consider a two-dimensional system (a thin elastic film) but the crack is 
in the present case two-dimensional and runs parallel to the plane of the film. This 
can be envisaged as a delamination problem: we are dealing with a crack running along 
a weak interface where the film is tethered to a rigid substrate. Failure is induced by 
shear tractions that are applied to the free surface of the film. The considered geometry 
is found in several physical applications, ranging from abrasive wear of coatings and 
shear-induced delamination of thin films to geophysical problems such as the failure of 
slopes by shear band formation jTj and the initiation of snow slab avalanches by failure 
of weak interfaces in snow stratifications [H [9]. From a theoretical point of view, a 
motivation for studying this system lies in the fact that the two-dimensional crack has 
additional degrees of freedom which, as we shall see, may lead to behavior which differs 
from that of one-dimensional cracks. 

2. Formulation of the model 

We consider an elastic film of thickness D tethered to a rigid substrate. The film is 
modeled as a regular network of identical coil springs of unit length and stiffness /. The 
(discrete) coordinates of the network nodes are denoted by X, Y , the network size is 
L X L, and the origin of the coordinate system is located in the centre of the system. 
The interface with the substrate is represented by breakable coil springs which connect 
each node to a rigid support as shown schematically in Figure [H We assume that the 
coil springs restrict displacements to be in the X direction only, thus we are dealing 
with a scalar displacement field U := Ux- The system is driven by an external shear 
force E acting in the X direction. This force is applied uniformly to all nodes of the 
network. 

The leaf springs break at the uniform critical extension Uq = 2. Structural 
randomness of the interface is modelled as quenched disorder in the form of random 
variations of the fracture energies Wp or, equivalently, of the leaf spring stiffnesses K = 
2Wf/Uq = W-f/2. These are assumed as independent identically distributed random 
variables with WeibuU distribution, P(Wf) = 1 — exp[— (IVf/Wq)'']- The parameter Wq 
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Figure 1. Schematic illustration of the model studied in this paper: A thin elastic 
film is connected by a cohesive interface to a rigid substrate; the film is loaded by 
shear tractions applied to its free surface. The film is modelled as a regular network of 
identical coil springs, while the interface is represented by breakable leaf springs which 
connect the nodes of this network to a rigid support. 



is chosen such that the mean value {W-p) = WqV{1 + 1//3) = 1 where T{x) denotes the 
Gamma function. Different degrees of randomness can be introduced by varying the 



exponent /5 and thus the coefficient of variation CV = \JT{1 + 2//3)/(r(l + — 1 

of the distribution. A crack is introduced by removing all leaf springs over the cracked 
area. We consider two crack geometries: (i) a crack of width 2A in the X direction that 
is located in the system centre and spans the system in Y direction, and (ii) a circular 
crack of radius A in the centre of the system. 

We note that a scalar model very similar to the present one was studied by Zapperi 
et. al. [To]. The main differences between both models concern the crack geometry and 
driving method: Zapperi et. al,. consider straight cracks that cut from one edge of the 
system halfway into the interface, and drive these by imposing a growing displacement 
on the cracked edge of the film. This geometry could be adapted to study size effects 
in systems with straight cracks, but it is not well suited for circular crack geometries. 

In the limit of slowly varying displacement fields, the force equilibrium condition 
for our model can be written in the following continuum approximation: 

where T is the local characteristics of the leaf springs representing the interface: 

T([/, X, F) = <( 2 ' - ' (2) 

0, U >2. 
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We note that Eq. ([T]) can be obtained from a long-wavelength approximation to the 
continuum equations describing stress equilibrium in a thin film under shear loads. 
The non-dimensional variables of the present model can thus be related to the physical 
and geometrical parameters of the thin film problem. The derivation has been given 
elsewhere [21 [H] and we refer the reader to those previous studies for making the link 
with 'real' physical systems. 



3. Displacement fields, critical crack lengths and critical stresses for linear 
and circular cracks 



For analysing the behavior of cracks in finite disordered systems, it is useful to first 
consider the same crack geometry in a system without disorder {Wy{X,Y) = Wp) and 
of infinite extension. We envisage two crack geometries, namely (i) a linear crack of 
extension 2A in the X direction which spans the system in the Y direction, and (ii) a 
circular crack of diameter 2A. 

We first consider a linear crack. A critical crack is characterized by the fact that 
the displacement on the crack line (in this case the lines X = ±A) has the critical value 
U = 2. With the additional requirement that the displacement gradient must vanish at 
infinity, we find that the solution of Eqs. ([1]) and ([2]) is given by 

^:{A' 
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IXI > A. 



The critical crack length follows from the requirement that the displacement must be 
differentiable at iXl = A: 



A 



-A* 



(4) 



where A* = ^21 /Wp is a characteristic length which governs the exponential decay of 
displacement and stress fields ahead of the crack tip. We note that the present relation 

S(v4 + A*) oc (5) 

represents a scaling that is different from Griffith's classical relation Sa/A oc ^/W^ for 
bulk cracks (which also holds for linear cracks in thin sheets as studied by Alava et. al. 
[6]). The different scaling reflects the fact that in case of a thin film, stress relaxation 
due to crack opening is confined to a volume that is limited by the film thickness (which 
is unity in the present non-dimensional representation). The characteristic length A* 
enters the expression ([5]) as an extension of the actual crack length A, i.e., it acts in 
a similar manner as the process zone size Ap in Bazant's generalization of Griffith's 
relation to bulk cracks in semi-brittle materials IT2l. 



A + Apoc 
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For a circular crack we transform Eq. ([T]) to cylindrical coordinates. With 
R = + r2 have 

Again, for a critical crack we require the displacement to assume the critical value U = 2 
on the crack line R = A. The solution is then 

^iA'-R') + 2, R<A, 

= ^ 2S ^ / _ 2S N Ko{^W^/{2I)R) ^ ^ ^ (7) 

W^F V VTf/ i^o(^VrF/(2J)A)' ^ ' 

where ii'n is the nth-order modified Bessel function of the second kind. The condition 
for matching the derivatives at R = A becomes 



/ 2S\ / Ki(JWf/(2I)A 

U'(A) = il^ = ( 2 - 4? ) JWf/2I ) ^ 



This is readily analyzed in the limit A ^ 1 where our underlying continuum 
approximation is expected to be valid. The first and second order modified Bessel 
functions have the same asymptotic behavior and their ratio rapidly converges to unity. 
Thus, we obtain the following condition relating stress and critical crack radius: 

This relation is (but for a factor of 2) identical to that for a linear crack. In both cases, 
plots of the inverse stress vs. the critical crack length are expected to produce a straight 
line with positive A-axis intercept, 1/S = + a where a = A* j ^2W-pI = l/Wp. The 
slope of this plot is expected to be s = 1/\/2WfI for linear and s = 1/(2v^2Wf/) for 
circular cracks. 



(9) 



4. Simulations of disordered systems 

In our simulations of disordered systems the coil spring constants I were set to unity 
and leaf spring constants were assigned randomly from a WeibuU distribution with 
exponent /3 and mean value (K) = 1/2 ((Wp) = 1). System sizes were varied between 
L = 10 and L = 160. To introduce a linear crack of width 2A we simply removed all 
leaf springs within a strip of this width in the centre of the system. Thus the crack 
widths investigated are integer numbers. Circular cracks were introduced by choosing 
a number B, removing all leaf springs for which -\- Y"^ < B^, and then evaluating 
the 'radius' of the crack a.s A = ^Jn/tt where N is the number of removed springs. For 
each set of parameters L, A) we performed - depending on system size - typically 
10^ simulations with different realizations of the quenched disorder. 

Simulations were performed in the standard manner for a random spring model. We 
solve the force equilibrium equations in the small-displacement approximation for unit 
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applied stress and identify the leaf spring with the largest displacement U . The external 
stress is then multiplied with S = 2/f/ such that this spring has the critical displacement 
U = 2, and the leaf spring is removed. We then again solve the equilibrium equations 
and identify the leaf spring with the largest displacement. If this exceeds the critical 
level, this spring is also removed and the process is repeated until a new equilibrium 
configuration is found where U < 2 for all sites where leaf springs are still in place. We 
then again identify the leaf spring with the largest displacement, multiply the external 
stress with 2/11, remove the spring, and so on until global failure of the interface occurs, 
defined in this case as the failure of all leaf springs. 

Results are shown in Figure [2] for circular cracks, and in Figure [3] for linear cracks. 
For circular cracks, two different types of behavior can be observed. At small crack 
radii, the failure stress is approximately independent of the crack radius A and has 
the same value as in a crack-free system. This failure stress has a logarithmic system 
size dependence as shown in the inset of Figure [31 At large crack radii, on the other 
hand, the failure stress is system size independent but depends on the crack radius. 
This dependence is similar to what we expect for a system without disorder: the plot 
of 1/E vs. A shows a linear dependence. The slope s and axis intercept a of this linear 
regime do not depend on system size. They exhibit only a weak disorder dependence as 
shown in the inset of Figure [2l The numerical values are close to the values a = 1 and 
s = l/{2\/2) expected in the continuum limit of a disorder-free system. 




1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 

critical crack radius A 



Figure 2. System failure stress versus crack length for circular cracks in systems of 
different size; (3 = CV = 1. Inset: Dependence of slope and axis intercept of the linear 
regime on the coefficient of variation of the fracture energy distribution. 

The failure stresses of systems with linear cracks shown in Figure |3] show different 
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behavior. Here, the two distinct regimes manifest from Figure [2] are absent. Instead, 
the plots of 1 /S vs. A are nearly linear for all crack sizes, and the logarithmic system 
size dependence of the failure stresses persists at all crack lengths. The slopes of the 
plots do not depend on system size. Their value is practically the same as for circular 
cracks, in good agreement with the deterministic result. 




crack length 2A 



Figure 3. System failure stress versus crack length, for linear cracks in systems of 
different size; /? — CV — 1. Inset: System size dependence of the strength of the 
crack-free system. 



5. Discussion and Conclusions 

The behavior we observe for circular cracks is analogous to the behavior of linear cracks 
in thin sheets described by Alava et. al. [6j. The cross-over between a system-size 
dependent but crack-size independent regime at small crack sizes, and a crack-size 
dependent but system-size independent regime at large crack sizes, can be understood 
in terms of a transition between crack-nucleation-controUed and crack-propagation- 
controlled failure. If the radius of the pre-existing crack is small, the stress needed 
for its propagation may be high enough to nucleate another critical crack in another 
part of a sufficiently large system. In this case, failure is not controlled by the initial 
crack but by crack nucleation in a 'weakest spot' elsewhere. This has been shown to 
give, in a thin-film geometry similar to the present one, rise to a logarithmic system 
size dependence of the failure stress [11]. If the radius of the pre-existing crack is large 
enough, on the other hand, failure is governed by propagation of this crack and does not 
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depend on system size. The crack length dependence of the failure stress matches the 
expectation from deterministic fracture mechanics: We find a linear dependence of the 
inverse failure stress on crack length, 1/S = a+sA. The parameters a and s of the linear 
fit exhibit a weak, approximately linear dependence on the coefficient of variation of the 
failure energy distribution. The slope s, which in deterministic fracture mechanics is 
proportional to s cx 1/ \/Wf, decreases slightly with increasing disorder. This indicates 
a disorder-induced increase of the effective fracture energy, which can be understood 
as a consequence of crack front roughening and crack pinning as discussed in |2]. The 
axis intercept a, by contrast, increases with increasing disorder. This behavior cannot 
be understood in terms of an increase in effective fracture energy. Rather, we may 
interpret the increase of a (or of A*) as a consequence of an advance of the crack beyond 
the initial crack line that is associated with crack front roughening under subcritical 
loads, an effect that becomes more pronounced as the disorder increases [2]. 

How can we understand the different size-dependent behavior of linear cracks? The 
obvious difference between the two geometries is that, for a crack spanning the entire 
system, the length of the crack line is equal to the system size whereas for a circular 
crack, it is proportional to the crack radius. As demonstrated in [2], crack lines in 
interface fracture exhibit self-affine roughening, but the regime of correlated roughening 
is confined to scales less than the crack width A. Above this scale, the critical crack can 
envisaged as a sequence of independently pinned segments of length ~ A. Supercritical 
crack propagation can be triggered by the depinning and expansion of any of these 
segments. Accordingly, the system strength is governed by the weakest link along the 
crack line. The validity of this weakest-link argument can be directly demonstrated 
by analyzing the probability distributions of failure stresses: If the argument is valid, 
the distribution P|"'-'(S) of the minimum failure stress determined from n independent 
simulations of systems of size L must be equal to the distribution PniiX) of the failure 

in) 

stress of a system of size nL. The n-th extremal distribution is calculated from 
the underlying distribution using the basic relation of extreme value statistics: 
P'i^\^) = !-[!- Pl(S)]". Results of a comparison between PnL and P^"'' for different 
values of L and n are shown in Figure |H It is evident that the weakest-link argument 
works quite well. We note, however, that the same is not true for system sizes that are 
comparable to the crack length A (such as L = 20 for the parameters used in Figure 
4). We attribute this to the correlated nature of the dynamics on scales < A, which 
precludes the use of simple weakest-link arguments. 

Establishing the nature of the associated size effect requires understanding of the 
functional nature of the failure stress distribution. Fitting the observed distributions 
indicates that, among the stable extreme value distributions, only the Gumbel 
distribution P(S) = 1 — exp[— exp(C(S — Eq)] does an acceptable job. Evidently, 
this implies a logarithmic size effect as observed in the simulations. 

In conclusion, the present investigation demonstrates that the additional degrees 
of freedom associated with 2D as compared to ID cracks may give rise to additional 
types of size dependent behavior. For 2D cracks as investigated in this study, the 
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Figure 4. Comparison of the probability distribution of failure stresses for a system 
of size L = 160 with the second extremal distribution for a system of size L = 80 and 
the fourth extremal distribution for a system of size L = 40. Parameters: crack length 
2 A — \Q, P — CV = 1. Full line: best fit using a Gumbel distribution with parameters 
C = 134.42, Eo = 0.2197. 



process zone width may be associated with crack front roughening rather than with 
preferential accumulation of damage ahead of the crack tip. The latter mechanism, 
which characterizes the process zone of linear cracks in 2D sheets where it leads to 
screening of crack-tip stress concentrations [6], is irrelevant in the present case where 
elastic couplings are short ranged: In the present geometry, stress and displacement 
fields ahead of the crack decay exponentially even in the absence of disorder, hence, 
the crack has negligible influence on its wider surroundings. In case of linear system- 
spanning cracks, failure is governed by the depinning of a crack-line segment with a 
width of the order of the crack length. Since system failure depends on the pinning 
strength of the weakest segment along the crack line, this gives rise to a size dependence 
of the failure stresses at all crack lengths. The nature of the failure stress distributions, 
as well as direct determination of the mean failure stresses for different system sizes, 
indicate a logarithmic size effect. 

Of course, a most interesting question concerns the size- dependent behavior of 2D 
cracks in bulk systems where both crack front roughening and damage accumulation 
ahead of the crack tip are likely to operate simultaneously. The interplay of these 
two mechanisms, the concomitant process zone characteristics, and the resulting size 
effects remain a task for further studies relying on much more substantial deployment 
of computer power. 
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